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A rich variety of order parameter manifolds of multicomponent Bose-Einstein conden- 
sates (BECs) admit various kinds of topological excitations, such as fractional vortices, 
monopoles, skyrmions, and knots. In this paper, we discuss two topological excitations 
in spinor BECs: non-Abelian vortices and knots. Unlike conventional vortices, non-Abelian 
vortices neither reconnect themselves nor pass through each other, but create a rung between 
them in a topologically stable manner. We discuss the collision dynamics of non-Abelian vor- 
tices in the cyclic phase of a spin-2 BEC. In the latter part, we show that a knot, which is a 
unique topological object characterized by a linking number or a Hopf invariant [773 (iS 2 ) = Z], 
can be created using a conventional quadrupole magnetic field in a cold atomic system. 

§1. Introduction 

Topological excitations exist in a wide variety of systems, such as liquid crystals, 
superfluids, superconductors, and our universe. Among them, Bose-Einstein conden- 
sates (BECs) of dilute atomic gases^® are the ones of the ideal testing grounds for 
investigating topological excitations. The novelty of this system lies in the unprece- 
dented controllability, which enables us to study static and dynamic properties of 
topological excitations. Various experimental parameters such as the temperature, 
atom number, trapping potential, and even the interatomic interaction can be con- 
trolled with high precision. 

When a BEC is trapped in a magnetic potential, the spin of each atom is oriented 
along the direction of the local magnetic field. The spin degrees of freedom are there- 
fore frozen and the BEC is described by a scalar order parameter. In such a system, 
the order parameter manifold is U(l), which can host quantized vortices. By using 
high-resolution imaging techniques, nonlinear dynamical phenomena such as vortex 
nucleation and vortex-lattice-formation dynamics,® and quantum turbulence® have 
been directly observed. 

On the other hand, when a BEC is confined in an optical trap, the direction of 
the spin can change dynamically due to the interparticle interaction. Consequently a 
BEC of spin-/ atoms is described with the (2/ + l)-component order parameter. In 
this case, the order parameter manifold is generally larger than U(l), and therefore, 
various nontrivial topological objects can be accommodated in a spinor BEC. 
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In this paper, we briefly overview the symmetry properties of the spinor BEC, 
and discuss two topological excitations: non-Abelian vortices (Sec. [3|) and knot soli- 
tons (Sec. Hi. 

§2. Symmetry property of spinor BECs 

We consider a system of spin-/ atoms under zero magnetic field. For simplicity, 
we consider a spatially uniform system. The Hamiltonian of this system is given by 

H= [ dr^£v$>)Wm(r) 

+ J drdr' £ ^'(ry^r^rO^O^r), (2-1) 

mnm'n' 

where ip m is the field operator for an atom in the magnetic sublevel m = 0, ±1, • • -±/, 
and M is the atomic mass. In dilute gases of neutral atoms, the main interatomic 
interaction is Van der Waals attraction and hard-core repulsion, which conserves the 
total spin T of two colliding atoms. Each scattering channel of total spin T can be 
described with the s-wave scattering length ajr, and V™71,(r,r') is then given by 

V™, n n ,(r,r') = 5(r-r>) £ ^-a T £ (fmfn\TM)(TM\fm'fn'), (2-2) 

.F=Q,2,- ,2/ M=-T 

where (fmfn\J-Ai) are the Clebsch-Gordan coefficients. The scattering amplitude 
for odd J- vanishes due to the Bose symmetrization. 

Since the interaction conserves the total spin and the number of particles, the 
Hamiltonian (|2-ip has the SO(3) rotational symmetry in spin space and the £7(1) 
global gauge symmetry, i.e., the full symmetry of this system is given by 

G = SO(3) F x £7(1) , (2-3) 

where subscripts F and (f> denote the spin and gauge, respectively. An element of 
G can be described with the combination of the Euler rotation e~ lFza e~ lFv ^e~ lFz,y 
and the gauge transformation e 1 ^ , where F n (/j, = x,y,z) is the /i component of the 
spin operator and a, f3, 7 are the Euler angles. While the system has the symmetry 
of G above the transition temperature, the symmetry G is broken to its subgroup 
H below the transition temperature. The subgroup H is referred to as the isotropy 
group. The order parameter manifold is defined as the coset space R = G/H. 

In the case of spinor BECs, several ground-state phases arise depending on 
the scattering lengths. In Table HI we show the representing order parameters of 
the ground state that appear in spin-1 and 2 systems, together with their isotropy 
groups H. The symmetry of the order parameter can be visualized by plotting the 
order parameter 

/ 

#(M)= E Yfm(e,^m, (24) 
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where Yf m (6, cp) is the rank-/ spherical harmonics with (0, ip) denoting the direction 
in the spin space. Figure Q] shows the surface plots of the ground-state order param- 
eters in spin-1 and spin-2 BECs. In the case of spin-1 BECP® the ferromagnetic 
order parameter has an SO (2) = U{\) symmetry axis as shown in Fig. [I] (a) (i). This 
corresponds to the fact that the order parameter ip F = (1,0,0) is invariant under a 
spin rotation e - *^ 2 ^ followed by a gauge transformation e 1 ^ , i.e., the ferromagnetic 
phase has the spin- gauge coupled U(l) symmetry. To clarify this particular feature 
of the coupling between spin and gauge degrees of freedom, we denote this symme- 
try as U(1)f z +c/>- On the other hand, the polar state shown in Fig. [1] (a) (ii) has 
Z2 symmetry in addition to the U{\) symmetry; ipp = (0,1,0) is invariant under 
e i-K e -iF x -K anc j g -iF z 0_ Therefore, the isotropy group H of the polar phase is given by 
the spin-gauge coupled dihedral group (Doo)f z ,0 = U(1)f z x (Z2)f,</>> where semidi- 
rect product x implies that when the nontrivial element of Z2 acts on an en element 
g € U(l), g changes to g . 

The symmetry for the spin-2 order parameters can be investigated in a similar 
manner, and the results are summarized in Table HI There are three mean-field 
ground states; ferromagnetic, cyclic and T)2™^i an( j quantum fluctuations divide 
the region of the D2 phase into uniaxial nematic (UN) and biaxial nematic (BN) 
phases.® As shown in Fig. [T] (b), the cyclic, D2, UN, and BN phases have discrete 
symmetries. In particular, the isotropy groups H for the cyclic, D2, and BN phases 
are non-Abelian. In these phases, vortices exhibit nontrivial collision dynamics as 
discussed in the next section. 



phase 
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MG/H) 


ir 3 (G/H) 
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(1,0,0) 


U{1)f z -4> 
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(Z)f 
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(T* x Z) FA 
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D2 


(a, 0, b, 0, 0) 


(D 2 )f z ^ 


(Q x Z) F , 





{Z)f 


UN 


(0,0,1,0,0) 




(Z 2 )f x (Z)* 


(Z)jp 


{Z)f 


BN 


(1,0,0,0,1) 




(Dt x Z) Ft4 , 





(Z) f 



Table I. Representing order parameter ip, isotropy group H of the order parameter, and the first, 
second and third homotopy groups of G/H for the ground states of spin-1 and 2 BECs. F, 
P, C, UN, and BN denote ferromagnetic, polar, cyclic, uniaxial nematic, and biaxial nematic, 
respectively. D n is the dihedral group, T is the tetrahedral group, and Q is the quaternion group. 
T* and D\ are defined by SO(3)/T 3* SU(2)/T* and SO{2,)/D A SU(2)/Dl, respectively. 

The n-dimensional topological objects are classified with the nth homotopy 
group ir n (G/H), which are also listed in Table HI Due to the rich variety of the or- 
der parameter manifolds, various kinds of topological excitations can exist in spinor 
BECsPH29 

§3. Non-Abelian vortices 

The topological charge of a vortex determines how the order parameter changes 
as one circumnavigates a loop surrounding the vortex. For example, if the phase of 



4 



Y. Kawaguchi, M. Kobayashi, M. Nitta, and M. Ueda 



(a) spin 1 (b) spin 2 




Fig. 1. Surface plots of the order parameter defined in Eq. (|2-4|l . where the gray scale on the 
surface represents argi?. F, P, C, UN, and BN denote ferromagnetic, polar, cyclic, uniaxial 
nematic, and biaxial nematic phases, respectively. The arrows labeled with (7(1) indicate the 
continuous symmetry axes, and those labeled with C n denote the n-fold symmetry axis. For 
example, in (a) (ii), the order parameter is axisymmetric around the U(l) axis, and has 2-fold 
symmetry about the C2 axis; when we rotate the order parameter about the C2 axis by ir, the 
phase of the order parameter changes by n, and to compensate this phase change, the gauge 
transformation e ln should be applied. 



the order parameter of a scaler BEC changes as e m< ^, the topological charge of this 
vortex is said to be n. In a similar manner, the vortices can be characterized with 
an operator which is acted on the order parameter as one goes around the vortex. 
To satisfy the singlevaluedness the order parameter, the operator should keep the 
order parameter invariant. 

We have shown in Ref. [18]) that the cyclic phase of a spin-2 BEC can host non- 
Abelian vortices. To understand this, we note that the order parameter of the cyclic 
phase is invariant under the following 12 elements of the tetrahedral group 1, 
I x = e iF ^, I y = e iF v\ l z = e iF *\ C = e ^i/3 e -2^HF x +F v +F z )/3^ > £ 2 , ^ ^ 
I Z C, I X C 2 , I y C 2 , and I Z C 2 , where we choose lobes in Fig. [1] (b) (ii) as x, y and 
z axes. These operators are regarded as topological charges of vortices in a cyclic 
BEC, and therefore, the vortices are non-Abelian. Furthermore, these operators are 
classified into four conjugacy classes: 

• (I) integer vortices: {1}; 

• (II) 1/2 - spin vortices: {I x , I y , I 2 }; 

• (III) 1/3 vortices: {C, I x O,I y C, hC}\_ 

• (IV) 2/3 vortices: {C 2 , I X C 2 , I y C 2 , I Z C 2 }. 

By acting an element of T, topological charges in the same conjugacy class trans- 
form into one another, which means that the topological charges are not uniquely 
determined. 

The non-Abelian characteristics of the vortices manifest themselves most dra- 
matically in the collision dynamics. In general, when two vortices collide, they 
reconnect themselves, pass through [Fig. E(d)], or form a rung that bridges the two 
vortices [Fig. [2(b), (c)]. When two Abelian vortices collide, all these three cases are 
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Fig. 2. Collision dynamics of two vortices, (a) Initial configuration, where A and B represent oper- 
ators that generate the corresponding vortex (a set of spin rotations and gauge transformations 
for the case of a cyclic BEC). The vortex on the bottom right, which is connected to B, is 
identified as ABA" 1 . The configuration in (a) is topologically equivalent to (b) and (c), where 
a rung is formed. If the vortices are a pair of vortex and anti- vortex, i.e., A = B" 1 , the rung 
in (b) disappears, giving rise to reconnection, whereas the rung in (c) corresponds to a doubly 
quantized vortex. If A and B are commutative, passing through is also possible because the 
configurations of (a) and (d) will then be topologically equivalent. However, when A and B are 
not commutative, the collision always results in the formation of a rung. 

possible, and one of them occurs depending on the kinematic parameters and ini- 
tial conditions. However, when two non-Abelian vortices collide, only a rung can 
be formed, whereas reconnection and passing through are topologically forbidden 
because the corresponding operators do not commute with each other. In fact, the 
nonvanishing commutator of the two operators gives the operator of the rung vor- 
tex. In Ref . [T8j) , we have numerically simulated the rung formation dynamics in the 
spin-2 cyclic phase. Figure [3] illustrates a typical rung formation. 




Fig. 3. Numerical simulation for the collision dynamics of non-Abelian vortices in a spin-2 cyclic 
BEC. Reprinted from Ref.lTH)). 

When the core of a 1/3 vortex in the cyclic phase is filled with the ferromagnetic 
state, which is possible for a certain parameter set, it is possible to observe such 
dynamics of vortex lines by using a phase-contrast imaging technique that can detect 
local magnetization. 

§4. Knot soliton 

Next we consider topological excitations in higher dimensions. The spin-1 polar 
phase and spin-2 UN phase can host 't Hooft-Polyakov monopoles and knots, while 
other phases in Table [I] support three-dimensional skyrmions. Among them, knots 
are distinguished from other topological excitations, such as vortices, monopoles, 
and skyrmions, in that knots are classified by a linking number while the others 
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are classified by a winding number. Knots are characterized by mappings from a 
three-dimensional sphere S 3 to S 2 . The same topological objects are discussed in 
high energy physics, where knots are suggested to exist as stable solitons in a three- 
dimensional classical field theory.^ 1 Here we explain that knots of spin textures can 
be created and observed in the polar phase of a spin-1 BECPS 

The S 3 domain is prepared by imposing the boundary condition that the order 
parameter takes on the same value in every direction at spatial infinity. Although 
the order-parameter manifold for the polar phase is complicated Rj> = (SO(3)f x 
U(l)^)/(D 00 )p j( f > = (Sp x U(l) ( f > )/('I<2)F,<f>, neither U(l) nor Z2 symmetry contributes 
to homotopy groups in spaces higher than one dimension. Therefore, when we discuss 
topological objects in higher dimensions, the order parameter can be described with 
a unit vector d € S 2 which corresponds to the direction of the "dumbbell" in Fig. [JJ 
(a) (ii). Consequently, we have 7r3(i2p) = ^(S" 2 ) = Z. The associated integer 
topological charge Q is known as the Hopf charge: 



Q = ^2 J ^ X CijkFijAk, (4-1) 

where Tij = diAj — djAi = d ■ (did x djd)W^ Note that the domain (r) is three- 
dimensional, while the target space (d) is two-dimensional. Consequently, the preim- 
age of a point on target S 2 constitutes a closed loop in S 3 . Furthermore, the Hopf 
charge is interpreted as the linking number of these loops: if the d field has Hopf 
charge Q, two loops corresponding to the preimages of any two distinct points on 
the target S 2 will be linked Q times [see Fig. HI (a)]. Figure 21(b) shows an example 
of the d field of a polar BEC with Hopf charge I ra 

Knots can be created by manipulating an external magnetic field. In the presence 
of an external magnetic field, the linear Zeeman effect causes the Larmor precession 
of d, while d tends to become parallel to the magnetic field because of the quadratic 
Zeeman effect. Suppose that we prepare an optically trapped BEC in the m = state 
[i.e., d = (0, 0, 1) T ] by applying a uniform magnetic field in the z direction. Then, we 
suddenly turn off the uniform field and switch on a quadrupole field. Because of the 
linear Zeeman effect, d starts rotating around the local magnetic field, and therefore 
the d field winds as a function of t, resulting in a formation of knots. Figure [5] shows 
the creation dynamics of knots in an optical trap subject to the quadrupole field, 
where the upper panels show the snapshots of the preimages of d = —z and d = x 
and the lower panels show cross sections of the density for m = — 1 components on 
the xy plane. The density pattern in m = — 1 components is the smoking gun of the 
knots; a double-ring pattern appears that corresponds to one knot. As the d field 
winds in the dynamics, the number of rings increases. This prediction can be tested 
by the Stern-Gerlach experiment. 



*' Strictly speaking, the configuration in Fig. 2] is an unknot, since the preimage of one point on 
S 2 forms a simple ring which is unknotted. 
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Fig. 4. (a) Preimages of two distinct points on S 2 forming a link, (b) Spin configuration of a knot 
with Hopf charge 1 in a polar BEC, where the arrows show the d field of the polar phase. The 
solid and dashed curves trace the point where d points to x and —z, respectively, forming a 
link. Reprinted from Ref. |20]). 



(a) 0.5 T L 



(b) 1.1 T L 

.d - (1, 0, 0) r 



d = (0, 0, -1) 




(c) 2.2 T L 



(d) 3.3 T L 





Fig. 5. Dynamics of the creation of knots in a spherical optical trap under a quadrupole magnetic 
field. Snapshots of the preimages of d = (0,0, — 1) T and d = (1,0, 0) T (top), as well as cross 
sections of the density for m = — 1 components on the xy plane (bottom). Reprinted from 
Ref. [H. 



§5. Concluding remarks 

We have discussed the symmetry properties of spinor Bose-Einstein condensates 
(BECs) and topological excitations in them. In particular, we have shown that 
when the BEC has the symmetry of non-Abelian group, the collision dynamics of 
the non-Abelian vortices are quite different from those of conventional vortices in 
a scalar BEC and superfluid 4 He. We have also shown that the spin-1 polar phase 
and spin-2 uniaxial nematic phase can accommodate a knot. The knot is unique in 
that topological excitations are classified with a linking number, while other topo- 
logical excitations are classified with a winding number. The superfluid helium-3 is 
often referred to as a testing ground to simulate our universe, because it can accom- 
modate various topological excitations similar to those known in gauge-field theory. 
The variety of topological excitations in spinor BECs is as rich as that in superfluid 
helium-3. The great advantage of spinor BECs is that direct observation and ma- 
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nipulation of the topological excitations are possible, while complicated analysis is 
required when using superfluid helium-3 if one is to identify a topological excitation 
from NMR signals. The controllability of spinor BECs would stimulate our imagina- 
tion, thereby promoting further developments in the study of topological excitations 
in the superfluid systems. 
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